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Abstract 

Here we report a theoretical model based on Green's functions and averaging techniques that gives ana- 
lytical estimates to the signal to noise ratio (SNR) near the first parametric instability zone in parametrically- 
driven oscillators in the presence of added ac drive and added thermal noise. The signal term is given by 
the response of the parametrically-driven oscillator to the added ac drive, while the noise term has two dif- 
ferent measures: one is dc and the other is ac. The dc measure of noise is given by a time-average of the 
statistically-averaged fluctuations of the position of the paramettic oscillator due to thermal noise. The ac 
measure of noise is given by the amphtude of the statistically-averaged fluctuations at the frequency of the 
parametric pump. We observe a strong dependence of the SNR on the phase between the external drive and 
the parametric pump, for some range of the phase there is a high SNR, while for other values of phase the 
SNR remains flat or decreases with increasing pump amplitude. Very good agreement between analytical 
estimates and numerical results is achieved. 

Keywords: parametric oscillations, stochastic differential equations, signal to noise ratio 



1 



I. INTRODUCTION 



Parametrically-driven systems and parametric resonance occur in many different physical sys- 
tems, ranging from Faraday waves [[D, inverted pendulum stabilization, stability of boats, balloons, 
and parachutes More recent applications in micro and nano systems include quadrupole ion 
guides and ion traps ['3l,opto-mechanical cavities [4], magnetic resonance force microscopy Q, 
tapping-mode force microscopy [|6|, axially-loaded microelectromechanical devices (MEMS) Q, 
and torsional MEMs [SJ, just to mention a few relevant applications. 

Parametric pumping has had many applications in the field of MEMs, which have been used 
primarily for measuring small forces and as ultrasensitive mass detectors since the mid 80's BH. 
An enhancement to the detection techniques in MEMs was developed in the early 90's that uses 
mechanical parametric amplification (before transduction) to improve the sensitivity of measure- 
ments. This amplification method works by driving the parametrically-driven resonator on the 
verge of parametric unstable zones. Rugar and Griitter [[TOll have shown ways, using this method, 
to obtain linear parametric gain. Furthermore, while they were looking for means of reducing 
noise and increasing precision in a detector for gravitational waves, they experimentally found 
classical thermomechanical quadrature squeezing, a phenomenon which is reminiscent of quan- 
tum squeezed states. Further experimental studies of parametric amplification appeared in [HTl . 
where a linear response around a limit cycle due to noise yields noise squeezing in a driven Duff- 
ing oscillator. Implementations of parametric oscillators in electronic circuits can be found in 
Refs. [[T2l[T3l . Parametric amplification started being studied in electronic systems in the late 50's 
and early 60's by P. K. Tien [|l4l, R. Landauer ||T5l, and Louisell [fT6l . It has been used for its de- 
sirable characteristics of high gain and low noise. Recent applications of parametric amplification 
in electronics can be found in Ref. [fTTl . 

The limits of parametric amplification due to thermomechanical noise on parametric sensing of 
small masses in nanomechanical oscillators have been studied in [[TSl . Although this work is quite 
broad the author does not provide estimates for the SNR as one approaches the first instability 
zone (e.g. by increasing the pump amplitude). Here, we study the effect of adding thermal noise 
to a parametrically-driven oscillator with the objective of studying the effectiveness of parametric 
amplification in the presence of noise. This one-degree of freedom model may be applied for 
instance to the fundamental mode of a doubly-clamped beam resonator that is axially loaded, in 
which case the one degree of freedom represents the amount of deflexion of the middle of the 

2 



beam from the equilibrium position. The present model can also be applied to the linear response 
of driven nonlinear oscillators to noise (such as transversally-loaded beam resonators), see for 
example Ref. [fTTI . One of us (A.A.B.) recently obtained analytical quantitative estimates of the 
amount of quadrature noise squeezing, heating or cooling in a parametrically-driven oscillator 
lfT9ll . We now use the Green's function approach, previously developed to solve the Langevin 
equation, aligned with averaging techniques, to obtain analytical estimates of the signal-to- noise 
ratio (SNR) in the parametrically-driven oscillator in the presence of both added noise and external 
sinusoidal drive. Here we show that for some values of phase (between the pumping and external 
drives), the signal grows faster than the fluctuations due to added noise, while for some other 
values of phase, the SNR is flat or decreases as the pumping amplitude grows, when one gets close 
to first instability zone of parametric resonance. 



II. THEORETICAL MODEL 



The equation for the parametrically-driven oscillator (in dimensionless format) is given by the 
damped Matthieu's equation 

X + oJqX = —•yx + Fp cos(2a;t) x, (1) 

in which 7 and Fp ~ 0(e), where e « 1. Since we want to apply the averaging method (AM) 
[|20l I2TII to situations in which we have detuning, it is convenient to rewrite Eq. ([!]) in a more 
appropriate form with the notation = to q — 10"^, where we also have ~ 0(e). With this 
substitution we obtain x + co'^x = —^Ix — 'jx + Fp cos(2a;t) x. We then rewrite this equation in 
the form x = y,y = —u?x + /(x, y, t), where f(x, y, t) = —ilx + Fp cos(2wt) x — 7?/. We now 
set the above equation in slowly-varying form with the transformation to a slowly-varying frame 
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After an application of the AM (in which, basically, we filter out oscillating terms at and near 2u} 
in the above equation), we obtain 
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where the functions U{t) and V{t) were replaced by their slowly- varying averages u(t) and v{t), 
respectively. The averaging theorem [21 J tells that these two sets of functions will be close to each 
other to order 0(e) during a time scale of 0(l/e) if they have initial conditions within an initial 
distance of 0(e). So by studying the simpler averaged system, one may obtain very accurate 
information about the corresponding more complex non- autonomous original system. With the 
transformations u{t) = e^^^^'^u{t) and v{t) = e~'^*/^f (t), we obtain 
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Upon integration of Eq. ([3]), one finds the solution 
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where k = — 5^, /3 = —Fp/Au, and 5 = VL/2uj. Hence, we find that the first parametric 
resonance, i.e. the boundary between the stable and unstable responses, is given by 



(5) 



This result is vahd for a; ~ even in the presence of added noise. In Fig. ([T]) we find very good 
agreement between results obtained from numerical integration of Eq. ([!]) and the boundary given 
by the averaging technique. 

We now will investigate the effect of added thermal noise on the parametric amplification mech- 
anism [[10l[nil22l|23l. We start by adding noise to Eq.Q and obtain 
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'jx + Fp cos{2ut) X + R(t), 



(6) 



where R{t) is arandom function that satisfies the statistical averages (-R(t)) = and {R{t)R{t')) = 
2T^5{t — t'), according to the fluctuation-dissipation theorem [|24ll . T is the temperature of the 



heat bath in which the oscillator (or resonator) is embedded. Once we integrate these equations 
of motion we can show how classical mechanical noise squeezing, heating, and cooling occur. 
We now summarize the method developed in [fT9l to analytically study the parametrically-driven 
oscillator with added noise, as given by Eq. ([6]). 

A. Green's function method 

The equation for the Green's function of the parametrically-driven oscillator is given by 

92 , d 



G{t,t') =6{t-t'). (7) 



Since we are interested in the stable zones of the parametric oscillator, for t < t' G{t,t') = and 
by integrating the above equation near t = t', we obtain the initial conditions when t = t' + 0"*", 
G{t, t') = and t') = 1.0. Using the Green's function we obtain the solution x{t) of Eq. Am 

in the presence of noise R{t) 

/oo 
dt' G{t,t')R{t'), (8) 
■00 

roc Q 

x{t) = v{t) = Vh{t) + j dt' -G{t, t')R{t'), (9) 

where Xhit) is the homogeneous solution, which in the stable zone decays exponentially with 
time; since we assume the pump has been turned on for a long time, Xh(t) = 0. By statistically 
averaging the fluctuations as a function of time we obtain 
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where t = t — t'. 

Although equation (|6]) may be solved exactly by using Floquet theory and Green's functions 
methods [1251 . one obtains very complex solutions. Instead, we find fairly simple analytical approx- 
imations to the Green's functions and, subsequently, to the statistical averages of fluctuations using 
the averaging method. We then use the solution of the system of coupled differential equations (|4]), 
where the initial conditions att = t' are given by u{t') = — sin (cut') /a; and v{t') = — cos{ujt')/u 



and obtain the approximate Green's function is 
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for t > t' and G{t, t') = for t < t' . In the stable zone of the parametrically-driven oscillator, 
when 1/3 1 > \5\ , we can rewrite the Green's function replacing the initial conditions and using 
simplifying trigonometrical identities. The change of variables t' = t — t leads to 
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We notice that by varying the pump amplitude Fp and the detuning f2, we can create a continuous 
family of classical thermo-mechanical squeezed states, generalizing the experimental results of 
Rugar and Griitter [fTOll . An estimate of the time average of the thermal fluctuations, when \(3\ > 
\5\,is given by 
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where the integrals are given by 
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A time-averaged estimate of the statistically averaged thermal fluctuations in velocity, when > 
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\S\,is given by 
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An estimate of the statistically averaged thermal fluctuations, when \(3\ > \6\,is given by 
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The remaining coefficients of eq. ( 16 1 are given by 
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Notice that when one gets close to the zone of instability we obtain a far simpler expression for 
the average fluctuations. It is given approximately by 
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III. LINEAR PARAMETRIC AMPLIFICATION 



In this section we study the parametric amplification of an added ac signal near the onset of the 
first instability zone of Matthieu's equation. The equation is given by 



X + U X 



-VLx — 7x + Fp cos(2ci;t) x + Fq cos(a;t + 



After doing averaging we obtain 
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The fixed points are given by 
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The gain of the amplifier is defined in [|26l as 
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where |X| = y/u*"^ + v*"^ and the pump off means Fp = 0. Rugar and Griitter [fTOl have studied this 
amplification process experimentally and also analytically via a perturbative method described by 
Louisell [fT6|. Although, their results agreed well with their experimental data, we believe that we 
can increase the applicability of mechanical parametric amplification of small signals by applying 
the averaging method and allowing for detuning. We also compare our analytical estimates of 



gain to the gain obtained from a full numerical integration of the equations of motion (18). The 



numerical gain is given by the expression in Eq. (21 1 with the analytical fixed-point values replaced 
by the numerical fixed points of the first-return Poincare map obtained from the integration of 



Eq. ( 18 1 after transients died out, i.e. x{nTy + y{nTY /u"^, where T = n /oj and n » 1. 
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IV. SIGNAL TO NOISE RATIO 



Following the previous definition of gain, we define a measure of the SNR as 

SNRo = 10 log 



(22) 



where the fixed points u* and v* are given by Eq. ([20]) and the time- averaged thermal fluctuations 



x^) are given by Eq. ( 14 1. Near the first instability zone we can write down this expression, with 

Fp [n cos(2</)) + 7w sin(2(/))] ' 



the help of Eq. ( 17 ), approximately as 

SNRo = 10 log 
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(23) 



Another measure of the SNR is given by comparing the signal intensity with the noise level at the 
same frequency, that is at 2u. This is given by the expression 



SNR 



'2u 



10 log 



U*2 + t;*2 



A 



2u}\ 



(24) 



2ui\ 



where the coefficients A2uj and B2U1 are defined in Eq. ( [T6] ). By dimensional analysis one notices 
that A2U) and i?2w have the same dimensional units as u*"^ and f *2. Near the first instability zone 
we obtain a simple estimate for this SNR measure, namely 
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V. RESULTS AND DISCUSSION 



In Fig. (|2]) we plot the Green's functions obtained directly from the numerical integration of 
Eq. ^ alongside analytical approximation results given by Eq. ( 12), if > or by Eq. (13) 



if 1/3 1 < \5\. We obtain excellent agreement between the two methods, what implies that our 
analytical estimates of {x'^{t)) are accurate. The numerical integration was performed using a 
RK4 algorithm with a time step given by vr/ (5 12a;). 

In Fig. (|3]) we obtain excellent fitting between analytical results for gain obtained by the aver- 
aging method in Eq. ( [2T| ) and numerical results given by the fixed point of the first-return Poincare 
map of Eq. ( [T8] ) obtained from the integration after transients died out. 

In Fig. (|4]) a comparison between numerical and analytical estimates of gain as a function of 
pump amplitude Fp is shown. Tow different values of phase are depicted. One observes a very 



strong dependance on phase between the pump and the external additive drive. One should obtain 
a divergence in gain as the boundary between stable and unstable response is reached. Again very 
good estimates are obtained. 

In Fig. ^ we show a logarithmic plot of the dc component of the mean square displacement 
over the heat bath temperature. The steep rise of the curve as the pump amplitude is increased 
indicates that the noise is also amplified by the parametric oscillator. 

In Fig. ([6]) we show logarithmic plots of the amplitude of the signal (amplitude of response of 
the parametric oscillator due to the external ac drive) over the noise (here the dc component of 



the mean square displacement). Note that the simple estimates given by Eq. (23) give very good 
approximation as one gets close to the first instability zone in parameter space. 

In Fig. (|7]) we show logarithmic plots of the amplitude of the signal (amplitude of response of 
the parametric oscillator due to the external ac drive) over the average fluctuations amplitude that 
oscillate at frequency 2uj (here the ac component of the mean square displacement). 
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VI. CONCLUSION 



The ability to reduce the influence of thermal noise in parametric amplifiers (which includes 
the linear regime of MEMS devices and optomechanical cavities), can greatly improve the accu- 
racy and precision of measuring small masses and weak forces. Here we extended the theoretical 
work related to the seminal experimental research by Rugar and Griitter [fTOl . For a long time 
since its publication there has not appear in the literature a sound explanation based on stochastic 
dynamics of the essential features of the classical thermal noise squeezing phenomena observed 
there. Recently, though, one of the authors (A.A.B) has proposed a stochastic dynamics model 
lfT9ll obtained by approximating the Green's functions of the parametric oscillator using averaging 
techniques to account for the observed experimental effects. Here we extend this work and give 
analytical approximate results for SNR in parametric amplification. Here we propose two different 
kinds of SNR (SNRq and SNR2^). These estimates of SNR give a measure of the effectiveness 
of the parametric amplifiers in the presence additive noise. The analytical results presented here 
confirm, that in both measures of SNR, that the parametric amplifier is indeed a very good ampli- 
fier, with sensitive amplification dependant on frequency, on phase and with low noise. Further 
refinement of our results may be achieved by including more details about the noise model, such as 
memory effects Il27ll28]| . by taking into account the coupling to a heat bath that could be made out 
of photons, as in radiation-pressure cooling, or via coupling to phonons. Further improvements 
of the accuracy of the predictions should be obtained by taking nonlinear terms into account, es- 
pecially when one gets close, in parameter space, to the first zone of instability. It is noteworthy 
to observe that our method may be applied to the linear response of nonlinear oscillators in the 
presence of both an ac drive and thermal noise. 

Finally, we note that this model can also be applied to the dynamics of ions in quadrupole ion 
guides or traps flSl or traps for neutral particles with magnetic dipole moments [|29ll . The presence 
of noise would indicate that the vacuum is not complete. One would obtain an estimate of the 
limits of mass spectroscopy in quadrupole ion guides, where the spectral limit is bounded by the 
ion guide length and the amount of noise present in the system. 
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FIG. 1 . (color online) Comparison between numerical and averaging method predictions for the boundary 
of the first instability zone of the damped parametrically-driven oscillator of Eq. ([T]). In the region above 
the green line lies the unstable zone obtained by numerical computation, while the region above the red 
line is the analytical prediction for the unstable zone. The numerical results are obtained by numerically 
calculating the corresponding Floquet multipliers, when at least one of them has modulus equal to 1. The 
averaging predictions are given by Eq. (|5]l. The fixed parameters of the equations of motion are 7 = 
0.1, cjQ = 1-0. These parameters were also used to obtain the results portrayed in the remaining figures. 
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FIG. 2. (color online) In the frames above we show several Green's functions with equally-spaced in 
time initial conditions in one given period of the parametric driving. They are vertically spaced for clarity, 
since all asymptotes are zero. In each frame we have a comparison between numerical results given by the 
numerical integration of Eq. (|7]l and the analytical approximate results given by Eqs. ( 12 1 or ( [T3] ). We have 
a.) Lo = 0.9 ujQ, b) = 1.0 ujq, and c) a; = 1.1 ujq- The initial values of the Green's functions are G{t, t') = 
and §^G{t, t') = 1.0 when t = t' + 0+. The pump amplitude used was Fp = 0.15. 
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FIG. 3. Comparison between numerical and analytical estimates of gain as a function of phase. The 



analytical gain is given by Eq. ( |2T] ). The numerical values are given by the expression in Eq. pT] ) with 
the analytical fixed-point values replaced by the numerical fixed points of the first-return Poincare map 
obtained from the integration of Eq. ( [TS] ) after transients died out and in accordance with the transformation 
in Eq. Q. Observe that both gain and absorption are reduced with detuning. 
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FIG. 4. Comparison between numerical and analytical estimates of gain as a function of pump amplitude 



Fp. The analytical gain is given by Eq. (21 1. The numerical values are given by the expression in Eq. (21 



with the analytical fixed-point values replaced by the numerical fixed points of the first-return Poincare map 



obtained from the integration of Eq. ( 1 8 1 after transients died out. 
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FIG. 5. (color online) Log plot of the do component of the mean square displacement, {x'^{t)), as given 
by Eqs. ([14]). In the linear parametric oscillator with thermal noise the temperature of the oscillator grows 
monotonically until it diverges at the boundary between stable and unstable zones. The approximating 



curves ai^e given by Eq. ( 17 1. 




FIG. 6. Plot of signal to noise ratio SNRq as defined in Eq. (22 1 along with corresponding approximations 
given by Eq. ( [23] ). 
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y=0.1, 03203=1.0, F2W=0.08 




FIG. 7. Plots of the signal to noise ratio SNR2a; as defined in Eq. (24 1, in which the signal is given by 
_j_y*2 ^jjj jj^g noise is given by the amplitude of squeezing at 2uj of the mean square displacement (x^ (t ) ) 



as given by \/p42^Jp~+~|^B^^, in which ^421.^ and i?2a; are defined in Eq. ( 16 1. Very precise approximations 



are also plotted, given by Eq. ( 25 I 
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